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ABSTRACT 

A conventional explanation of the dipole anisotropy of the cosmic microwave back- 
ground (CMB) radiation is in terms of the Doppler effect: our galaxy is moving with respect 
to CMB frame with ~ 600 km s~ x . However, as the deep redshift surveys fail to reveal a 
convergence of the large scale flow to zero at distances as large as d ~ H _1 15, 000 km s _1 
(Lauer & Postman, 1994), the uniqueness of the conventional interpretation has to be in- 
vestigated. A possible alternative might be a cosmological entropy gradient, as suggested 
by Paczynski & Piran (1990). We find that contrary to that suggestion a quadrupole 
anisotropy is generically of the same order of magnitude as the dipole anisotropy (or 
larger) not only for adiabatic but also for iso-curvature initial perturbations. Hence, the 
observed dipole cannot be explained with a very large scale perturbation which was initially 
iso-curvature. 



Subject headings: cosmic background radiation - cosmology - gravitation 
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1. INTRODUCTION 

The dipole moment of the CMB is usually interpreted to be the result of a Doppler 
effect caused by our motion with respect to CMB frame (cf. Kogut et al. 1993 and 
references therein). According to this interpretation the CMB, as well as galaxies (when 
averaged over a large enough volume) define the local standard of rest. Our galaxy, together 
with its neighbours moves with respect to the local standard of rest. When the velocities 
of galaxies in a large enough volume are measured they should be found to be at rest, 
naturally after the allowance is made for the overall Hubble expansion. 

The observations available so far fail to provide a clear support for this picture. The 
recent most troublesome result is that of Lauer & Postman (1994) who find that the 
frame defined with the 119 Abell clusters of galaxies within 15,000 km s _1 is moving at 
~ 700 km s _1 with respect to the local standard of rest as defined by the CMB. This 
scale is so large that it is difficult to accomodate within most currently available models 
of the formation of large scale structure in the universe (Strauss et al. 1994). This trouble 
persisted for many years (cf. Paczyhski & Piran 1990, hereafter PP, and references therein) 
and it justifies a search for alternative interpretations of the CMB dipole anisotropy. PP 
proposed that entropy gradient on a scale larger than the current horizon could give rise 
to a dipole moment while keeping the quadrupole unmeasureably small. 

The purpose of this paper is to demonstrate that the PP proposal was incorrect. PP 
used Tolman-Bondi cosmological model with no pressure to demonstrate that the very long 
wavelength density perturbations give rise mostly to a quadrupole anisotropy, while entropy 
perturbations show up as a dipole. The reason for ignoring pressure was the currently very 
small value of density and pressure due to CMB. The CMB was dynamically important 
only at redshifts larger than z eq ~ 10 4 , and it seemed safe to ignore it, or at least it seemed 
to be a reasonable first approximation. It turns out that was a conceptual mistake. The 
purpose of this paper is to demonstrate that no matter how small is the current contribution 
of radiation to the closure density any large scale iso-curvature (entropy) perturbation 
generically gives rise to a quadrupole CMB anisotropy which is larger than the dipole 
anisotropy. We demonstrate this in section 2 for a plane wave, and in section 3 for a 
modified Tolman-Bondi model. Finally, a simple qualitative demonstration of the fact 
that the initially iso-curvature (entropy) perturbation generically gives rise to a density 
perturbation (cf. Peebles, 1993, hereafter P93) is presented in section 4, together with a 
discussion of other possibilities. 
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2. PLANE WAVE PERTURBATIONS IN A FLAT UNIVERSE 

In the synchronous gauge (P93) the perturbed metric of the flat universe model has the 
form: 

ds 2 = c 2 dt 2 - a 2 (t)(5 jm - h jm )dx J dx m (2.1) 

where t is the cosmic time, a(t) the scale factor in the unperturbed model, Sj m is the 
Kronecker delta symbol and hjm(t,x l ) are the small perturbations to the metric. The 
spatial coordinates are Cartesian, Latin indices go through 1,2,3. We are interested in the 
growing, scalar modes of perturbations to the metric (Lifschitz, 1946). In the case of the 
plane wave perturbation with the wave vector k along the x 3 axis, the only nonvanishing 
components of hj m are the diagonal terms and hu = h 2 2 because of the symmetry. Thus 
there are two independent variables defining tensor hj m and we choose h = hu + /122 + ^33 
and H = hu — /133 for the purpose. 

In our model the matter consists of two independent components, which interact 
which each other by gravity only. The first component is non relativistic (NR) matter 
(for example baryons) and it is characterised by the present dimensionless density Qnr = 
Pnr/Pc, where p c = 3Hq/8tvG is the critical cosmological density, Hq is the Hubble's 
constant and G is the constant of gravity. The second component is ultra-relativistic 
(UR) matter (for example electromagnetic radiation) and it is characterised by the present 
dimensionless density Qjjr — Pur/ Pc- Since we are using the flat unperturbed cosmological 
model, we add the cosmological constant to be in agreement with unperturbed evolutionary 
equations: 

d\ 87tG , .1,9 . . 

a) =^(^ + ^) + 3 Ac (2 ' 2) 

< ^ = -^{t NR + t UR + ?,P UR ) + -Kc 2 (2.3) 

where dots mean the time differentiation, a = ptc 2 depict the full energy densities of 
various components and P$ - their pressures. We introduce dimensionless cosmological 
constant A = ^Ac 2 /Hq. In a flat model we have: 

tt NR + n UR + X = l (2.4) 

In some cases A may be negative. 

Throughout this paper we assume that all perturbations of all non-relativistic matter 
follow strictly the perturbations of baryons. We also assume that all perturbations of 
all ultra-relativistic matter follow strictly the perturbations of electromagnetic radiation. 
The ultra-relativistic and non-relativistic components are coupled to each other prior to 
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recombination, and they become decoupled when the universe becomes transparent to 
radiation at z rec ^> 1. 

The energy density fluctuations of different components are described by: 

S -£SSL = iNR S JS^ = U oa (2.5) 

Pnr Pur <j 

where Snr, Sjjr are the relative perturbations in particle density of the two componets. 
The perturbations in specific entropy 5s = 5S/S, where S is the entropy of ultrarelativistic 
fluid per particle of nonrelativistic fluid can be exppressed as 5 s = 5jjr — 5nr- (We neglect 
the entropy of the nonrelativistic fluid) . The velocity perturbation of each fluid has the 
form u a = (1/c, i^/ac) where v is the physical velocity measured by the synchronous 
observers . In the scalar mode of perturbations only the correction to the expansion, O, 
enters the equations: 

< = 3- + e e = -x (2.6) 

a a 

The dimensionless variable x is more convenient than O (Press & Vishniac, 1980, hereafter 
PV80) and we need Xnr and xur to characterize both components. 

Our set of equations describing the evolution of small perturbations to the metric 
and fluid variables is based on P93 and PV80. We use a new time coordinate (PV80) 
rj = ln(a/ao) = —ln(l + z), where ao is the present characteristic scale in the Universe 
and z is the redshift. The spatial gradients are already omitted in the equations, so they 
are valid for very large scale perturbations only. The equations are written for two fluid 
interacting only gravitationally (but see below): 

... faa A ,, 8nGa 2 ( . 8 „ \ , n _. 

h + ( ^2 + 1 ) h = —2~ [PnrSnr + -PurSur) (2.7) 



an 



8nGa 2 ( , 4 



H" + I — + 2 1 H' = ^- Ipnr5nr + -PurSur) - h' (2.8) 

&NR=(\h , -XNR) (2-9) 

&UR = (\ti-XuR) (2-10) 

. / da \ , 

Xnr + ( -rp + 1 J Xnr = (2.11) 

, / aa\ , 

Xur + ^tj J Xur = (2.12) 

where we have neglected sound velocity in the nonrelativistic component and for the rela- 
tivistic component we have already substituted 1/3 for its pressure to energy density ratio 
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and for its sound velocity square. The time derivatives of a can be substituted from Eqs. 
(2.2), (2.3). 

We solve our equations for different values of the universe model parameters varying 
the density of non-relativistic component in the limits 0.01 < Qnr < 1 and the density of 
ultra- relativistic component is in the range 10 -5 < fljjR < 1- The value of cosmological 
constant results from Eq. (2.4). 

We always set initial conditions long before the time of recombination and long before 
the time when the energy density of the ultra-relativistic component is equal to the energy 
density of the non-relativistic component, whichever comes first. In this early time the 
radiation dominated plasma behaves like a relativistic fluid and the unperturbed model 
evolves like a model with relativistic equation of state. (Various terms in the equations have 
different redshift dependence. While pur ~ (1+z) 4 and pnr ~ (1+z) 3 the terms including 
the cosmological constant do not depend on the redshift. Thus the effects of A ^ can 
only show when z is small.) In the early times the coefficients in the above equation set 
remain constant so one may find independent modes of perturbations (Lifshitz, 1946, P93 
and references therein). We consider two types of initial conditions. The first, adiabatic 
perturbation, has three possible modes in the relativistic Universe. The most natural of 
them , with growing rate proportional to a 2 , remains regular at t — ► and has vanishing 
velocity (xnr = 0, xur = 0). Another growing mode in the early epoch with amplitudes 
~ a 1 has nonvanishing velocity and irregular some of the metric components when t — ► 
(P93). Since the rate of instability in this mode is slower as compared to the previous 
one, it can not appreciably influence the present Universe unless the initial conditions are 
fine tuned. The decaying mode is of no interest: first it cannot produce any significant 
perturbation to the present Universe, second it is unphysical (PV80, Bardeen, 1980). For 
adiabatic perturbations we have o~nr = Sjjr initially and this condition is preserved if 
velocities vanish. (This mimics the coupling between ordinary matter and radiation, if 
required, without explicitly puting it into equations). 

The second, isocurvature type of perturbations (Peebles, 1987) is impossible in a single 
fluid model. With more than one fluid we are able to introduce a spatial dependence of the 
chemical composition of the matter not introducing perturbations to the energy density or 
to the metric. We just perturb the density of non-relativistic component, not changing the 
density of ultra-relativistic component. As long as the relativistic component dominates, 
no energy density perturbation arises and perturbations change slowly (5^^ ~ const, h', 
H', Sjjr << Snr)- When nonrelativistic component becomes dynamically important the 
growing energy density perturbation results. Examining Eqs. (2.9) and (2.10) we see, that 
the entropy perturbation 5s = Sur — 8nr remains constant. 

The two kinds of perturbations considered behave differently at the begining, but 
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become similar when the model becomes nonrelativistic. We are interested in the fluctua- 
tions in the microwave background caused by the perturbations. They arise between the 
surface of last scattering and the present epoch. To compare results in different models we 
normalize all gravitational instability calculations in such a way that the present amplitude 
of the density perturbations of non-relativistic component is the same and has the value 
5. According to the above remarks it is not surpraising that the influence of both kinds of 
perturbations on the microwave background is similar. 

To find the temperature of the CMB radiation in any particular direction on the sky 
n, we have to follow rays back in this direction. The coordinate distance travelled by a 
particle moving with the velocity of light is given by 

T («)= rW ,))=/'£*=r£*» (2.i3) 

J a J a 

The position of a photon, which is now at xq, coming from the direction n, was at the 
"time" r\: 

x^ = x + n(r - t v ) (2.14) 

where we have substituted rj = for the present epoch. To measures the present coordinate 
distance to the horizon. 

The temperature fluctuation of the CMB radiation coming to the observer at xq from 
the direction n on the sky is given as (Sachs and Wolfe, 1967): 

(2.15) 

where the subscript "rec" denotes the quantities measured at the epoch of recombination, 
at the place from which rays come. The first term coresponds to the temperature fluc- 
tuations in the plasma at the last scattering surface. It is equal to \5ur in the place 
of emission. For a plane wave with a wave vector k we define the directional cosine 
(j, = cos 9 = nk/k. In our case k is along the x 3 axis and phase factor along the ray 
changes like exp(ikx 3 + ikfJ,(ro — r)). k(ro — r rec ) ~ kro << 1 for superhorizon perturba- 
tions. Expanding and taking real part we have to the lowest order: 

XT' \ 1 

— ) =-- kr sin(kx 3 ) S UR (r] rec ) P\{p) (2.16) 

Pi(p) — jx is the Legendre polynomial of the first order and the above expression con- 
tributes to the dipole anisotropy of the CMB (this contribution we shall denote D Tec ). 
The quadrupole part (Q re c) is of still higher order being proportional to (kro) 2 . 

The second term is due to the Doppler shift caused by the difference in velocity of 
matter between recombination and the present epoch. Since we do not consider velocity 
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perturbations, v rec vanishes automatically. The observer peculiar velocity v b s may arise 
from the small scale perturbations to the gravitational field but is of no interest here. 

The third term is the result of gravitational field acting on photons (the Sachs - Wolfe 
effect). With our definitions of h, H we have: 

\ h' jm n>n m = l -ti Po(ft) - l -H' P 2 ( M ) (2.17) 

where P n (p-) are the Legendre polynomials .* 

The following of rays is required in the integration of Eq.(2.16). The gravitationally 
induced temperature fluctuations are given as: 

(^) 0*) = -^***" / {h'{r } )PM-2H'{r } )P 2 { l i)) e ik ^-r v ) dr] {2A8) 

V ± / SW " J Vrec 

For a superhorizon perturbation the exponent under the integral can be expanded. We 
get the series in the small quantity kro with coefficients being the products of Po and P 2 
with powers of \i. We limit ourselves to the dipole and quadrupole terms: 

#sw = -kr sm(kx s ) f (\h\ V ) - ^H'(r,)\ (l - -^ dr, (2.19) 

f° 1 

Qsw = - cos(kx s ) / -H'(?])dr] (2.20) 

One can define the density difference accross the horizon A = |toV5| ~ kr^d. For 
various contributions to anisotropy we have: 

-Dree = d Tec A Dg W = dsW A (2-21) 

and 

A 2 

Qrec = Qrec ~ j~ QsW = ?SW $ (2.22) 

In Fig.l we show the ratio of the dipole to quadrupole CMB anisotropy (D/Q px D^w/Qsw 
measuring it by the ratio of the density difference through the horizon to the density 
perturbation amplitude (A/5) which is a small quantity being the ratio of the present 
horizon size to the present scale of perturbation. The dipole anisotropy caused by the 
superhorizon perturbation is always smaller than the resulting quadrupole as can be seen 
on the plots. 



* In general case also the amplitudes of vector metric perturbations coupled to the 
spherical harmonics Y±\ and the tensor amplitudes times Y"± 2 would appear in eq.(2.17) 



The expansion of the perturbed model is not isotropic and the "Hubble constant" 
depends on the direction of measurement. Comparing the proper distance to a close 
object in the direction n with its velocity due to the expansion one gets: 



H(n) = -- h lm n l n m = - 
a 2 a 



1 - -U'(0, x 3 )P (p) + V(0, x 3 ) P 2 (p) 

D 3 



(2.23) 



where the variables are calculated at the present time. The monopole part is of no interest. 
Perturbations in the metric introduce the quadrupole anisotropy to the Hubble law with 
the relative amplitude: 



Q H = -cos{kx 3 ) H'(0) 



(2.24) 



As one can see, the Hubble anisotropy has the same phase of spatial dependence as the 
quadrupole anisotropy of CMB, but the amplitudes are defined by different quantities. 



3. SPHERICAL SOLUTIONS 

We follow here the approach to the spherically symmetric world models outsketched in 
the Appendix B of PP. As in the previous chapter and generally, when the inhomogeneities 
have the scale much larger than the horizon scale at the epoch of interest, we are going 
to neglect the influence of pressure gradients on the evolution of the model. We use the 
Bondi-Tolman metric: 



ds z 



c 2 dt 2 - X 2 (t, r)dr 2 - R 2 (t, r)d& 



(3.1) 



where r is a radial coordinate and dfl 2 is the line element on the sphere. (In fact the 
presence of pressure gradients forces one to use a more general spherically symmetric metric 
with g tt = A 2 {t, r), since the field equations say that the pressure gradient is proportional 
to the gradient of A(t,r), as shown by May & White (1967). Thus, strictly speaking, 
the synchronous, comoving coordinate system is impossible in the presence of pressure 
gradients.) Neglecting pressure gradients from the begining, we adopt line element (3.1) 
in our calculations. 

In the case of a single fluid one can choose the coordinate system (3.1) to be comoving 
with the matter, so the velocities vanish automatically. In the case of many fluids inter- 
acting only gravitationally it is possible that pressure gradients, not necessarilly equal in 
different components, may cause the relative motion of the fluids. But we are neglecting 
pressure gradients from the begining, so it is fair to assume that relative motions of the 
fluids are impossible. Thus it is possible to make the coordinate system comoving with 
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the matter. In that case the mixed components of the energy-momentum tensor vanish 
automatically (T tr = 0) and, as a consequence of the field equations one has: 

where W(r) is a free function (Bondi, 1947, PP). Equation (3.2) implies, that one can 
write down the field equation for R(t, r) in the form: 

tf = WHr)-l+ 2 -^$^ (3.3) 

R 

where 

nr fT 

m(t, r) = 4n pR 2 R r dr = 4n pR 2 XWdr (3.4) 

Jo Jo 

The quantity defined above is the gravitational mass in the sphere inside r. The second 

integral shows that it is the total matter density integrated over the proper volume with 

some "weighting" function W(r). Equations with the form of Eqs (3.3)and (3.4) are valid 

in the general case of spherically symmetric configuration with pressure (May & White, 

1967). Equation (3.3) has the form of energy equation for a test particle in the field of 

[possibly variable] mass m. While the density distribution inside the configuration defines 

the time dependent potential, the function W(r) sets the initial velocity of the particle. 

If the pressure vanishes, mass inside any radius r is preserved and Eq (3.3) becomes an 
ordinary differential equation with coordinate r playing the role of a parameter (PP) . In the 
general case of matter having an admixture of relativistic component (at least photons are 
such a component) one may solve the evolutionary equation for R(t, r) using the following 
approach. First we divide the configuration into a number of concentric, thin shells. The 
innermost spherical region behaves like a part of uniform solution. For this region one has: 

m NR(t) = mlfniUmt) m\j R {t) = m\j R {t init ) *"'*'? (3.5) 

n{t, n) 

where m]^ R , m]j R are the masses of all non- relativistic and ultra-relativistic compo- 
nents, respectively, inside the first zone. The same behavior is true for any other non- 
relativistic/ultra-relativistic components. During adiabatic expansion a relativistic fluid 
changes energy density as Pur ~ V~ 4 ' 3 where V is the proper volume of the region. Using 
Eq. (3.4) we have for a thin shell of matter between rj and r J+ i: 



J+ 1 (+\ ™J 



J+ 1 (+ \ ™J 



^UR V / '^UR\ ) — ^UR \pifiit) fT^u R \J , init, 



■rt [J'initi Tj-\-l) -*£ Kpiniti Tj ) 

R^t,r j+1 )-R^t, rj ) 



(3.6) 



where we implicitly assumed that W(r) can be treated as constant through a single zone. 
Using Eq. (3.5) one can obtain solution for R(t, r\) and m(t, r\). Knowing R(t,rj) and 
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m(t,rj) one can find m(t, 7j+i) as a function of R(t, rj+±) with the help of Eq. (3.6). 
Substituting this dependence sinto Eq. (3.5) one finds solution for R(t,rj + i). Recurently 
one can find the metric functions values and their derivatives on a grid. In other points 
metric can be found by interpolation. 

We start calculations early, when ultra-relativistic component of matter dominates 
the dynamics (pur » Pnr- At this stage we always set relativistic fluid density to be 
constant in space. The perturbation is in the free function W(r), which we borrow from 
PP: 

2 2 2 

W 2 M = !-£-£:, (3.7) 

'0 "+" ' '0 

where r$ » 1 sets the spatial size of perturbation. In our convention coordinate distance 
r ^ 1 corresponds to the present horizon size. The density of non-relativistic component 
is also set to constant in the case of adiabatic perturbations. For entropy perturbations 
we use the following shape of initial density of the non-relativistic component: 



Pnr ~ &nr 
Accordingly the specific entropy behaves like: 



r \ 



r < 2r (3. 



S(r) ~ P NR -—3 (3-9) 

We put observers at different places in the world models described above. The co- 
moving volume is defined as: 

mr)= «£l (3.10) 

Observer at any location (£ bs, r), where t bs is to represent the epoch of observation, can 
define the recombination epoch time t rec using the condition: 

V{t ohs , r) = (1 + ^ec) 3 ^(t rec , r) (3.11) 

where z rec ~ 10 3 is the redshift factor of recombination. Suppose we follow a ray, which 
goes along a null geodesies from the last scaterring surface at (t rec , r rec ) to the observer 
at (t bs7^obs)- The starting point of the photon depends on the direction n from which 
it arrives, r rec = r rec (n). The locally measured energies of a photon would be E rec and 
£^obs respectively, the co moving volumes at the locations - V rec and V^bs- Since initially 
the radiation temperature was constant through the space and it changes accordingly to 
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the law T ~ I/ -1 ' 3 , we can define it at any time and location. Taking into account the 
redshift of photons between last scattering surface and the observer we have: 

i 
T(n) = T local §^ (^») ' (3.12) 

where Ti OC ai is the locally measured, average temperature of the CMB. We compare results 
of such CMB temperature measurement in the direction radially in (Ti), radially out (T3) 
and in transverse direction (T2). The average temperature on the sky in our approximation 
is given as: 

Ti OC ai = ^(T 1 + 2T 2 + T 3 ) (3.13) 

We define the dipole and quadrupole anisotropy of the CMB as: 

D = yf^ Q= r 1 -2r 2 + r 3 (3l4) 

" J- local ^ -L local 

According to our calculations the dipole is always about 2 orders of magnitude weaker as 
compared to the quadrupole. The exception are the places, where the quadrupole vanishes 
locally, but such places are rare and the a'priori probability of making observations from 
there is low. 



4. DISCUSSION 



In sections 2 and 3 we presented a formal demonstration that a perturbation which is 
initially iso-curvature (i.e. the entropy is perturbed) while the ultra-relativistic component 
dominates, grows into a curvature (density) perturbation when the overall expansion of the 
universe makes the non-relativistic component dominant. The dynamics of the universe is 
believed to be currently dominated by the non-relativistic component (e.g. baryonic, or 
cold dark matter), while it was dominated in the past by the ultra-relativistic component 
(e.g. radiation, pairs, etc.). We can envision the expansion history of two universes 
which differ in the initial entropy per baryon, or any other measure of the ratio of ultra- 
relativistic to non-relativistic components, both universes being exactly flat, isotropic and 
homogeneous. * 

The two expansion histories are synchronized at to = 0. Initially, they are identical, 
with the scale factor increasing as t 1 ' 2 while the dynamics is dominated by the ultra- 
relativistic component. The expansion rate changes to a ~ t 2 ' 3 at the time t = t eq , 



* Similar reasoning can be found in Tolman, 1934 and P93. 
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when the non-relativistic component becomes dominant. This time is different for the two 
universes with the different ratio of the two components, and therefore, the late expansion 
rate is different too. This means that the scale factors for the two universes are not the 
same at any time t > t eq . If these two universes are just large parts of the same universe 
and they come into contact at t > t eq , the mis-match of their scale factors will create 
space curvature, which will give rise to a quadrupole anisotropy of the CMB. On Fig. 2 
we show the expansion histories of different parts of our spherical model. At the begining 
R(t, r)/r ~ t 1 ' 2 independent of position. In later times the differences become visible. 

Our conclusion is that the dipole anisotropy of the CMB cannot be explained with 
a very large scale entropy gradients in the universe, as proposed by Paczyhski & Piran 
(1990). If the dipole is cosmological in origin something else is needed. We see no 'natural' 
explanation. A trivial and entirely ad hoc and artificial 'explanation' can be offered. 
Imagine there are many ultra-relativistic components, with one having a very large scale 
fluctuation which is exactly in the opposite phase than radiation. There would be no 
overall change in the ratio of ultra-relativistic to non-relativistic componet, and no effect 
on the dynamics of the universe, as any effect of the radiation perturbation would be 
exactly balanced by the opposite effect due to another ultra-relativistic component, yet 
there would be a dipole in the CMB. There is no justification for such a proposal, and 
we mention it as an example of what may have to be considered if the Lauer & Postman 
(1994) result is confirmed and even extended to ever larger scales. 

This project was supported with the NSF grant AST 93-13620, NASA grant NAG5- 
2759 and KBN grant 2 P304 006 06. 
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FIGURE CAPTIONS 

Fig.l. The ratio of the dipole D to quadrupole Q anisotropy of the CMB as a function of 
dimensionless matter density Qnr- The ratio of the density difference across the horizon 
A to the density perturbation amplitude 5, which is a small number for superhorizon 
perturbations, serves as a unit of D/Q. (See text for explanations of these parameters). 
The curves are for Q UR = 10~ 5 , 10~ 4 , 10~ 3 , 10~ 2 , 10 _1 and 1 from bottom to top. The 
cosmological constant is given by A = 1 — Qnr — ^ur m each case and may be negative. 
(a) results for the adiabatic and (b) isocurvature perturbations. (See the text for details). 



Fig. 2. The evolution of expansion factor R(t, r)/r with time for r = 0.5, 1.0, 1.5 and 2.0 xro 
(bottom to top). We use the results obtained for our spherical model with Qnr = 0.1 and 
Qjjr — 10 -4 . That illustrates the different expansion rate in different parts of space. 



